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Abstract. - Coherent backscattering is due to constructive interferences of reciprocal paths 
and leads to an enhancement of the intensity of a multiply scattered field near its source. To 
observe this enhancement an array of receivers is conventionally placed close to the source. Our 
approach here is different. In a first experiment, we recover the coherent backscattering effect 
(CBE) within an array of sources and a distant receiver using time correlation of diffuse fields. 
The enhancement cone has an excellent spatial resolution. The dynamics of the enhancement 
factor is studied in a second experiment using correlation of thermal phonons at the same 
ultrasonic frequencies, without any active source. 



Coherent Backscattering Enhancement (CBE) is a consequence of interferences in multiple 
scattering or multiple reverberations of waves in a complex environment. It was first observed 
for electrons [1], then in optics [2,3] where it exhibited as a doubling of the backscattered 
energy in the incident direction. Later, it was applied to various fields of wave physics, like 
ultrasound [4-6], laser in cold atoms [7], and more recently to seismology [8]. In these fields, 
CBE was shown to be an accurate way to measure transport mean-free paths of the medium, 
a quantity that characterizes its degree of heterogeneity. The origin of CBE is the constructive 
interference between long reciprocal paths in wave scattering or reverberation. This enhances 
the probability for the wave to return to the source by a factor of exactly two in open media, 
and three in closed cavities, which results in the local energy density enhancement by the 
same factor [9,10]. In a reverberant body, it was proved that the CBE factor increases from 
two to three with a dynamics driven by the Heisenberg time Th (also referred to as the break 
time of the cavity) [6]. When the source and the receivers are placed inside the diffusive 
medium, the spatial extent of the cone is the order of the wavelength [6, 11]. In all of these 
active experiments, the observation of CBE requires that a receiver be placed close to the 
source. Unfortunately, the source-receiver configuration is sometime disadvantaging. First, it 
may break the reciprocity condition if the source and the receiver have not exactly the same 
orientation [12]. Second, placing a source close to the receiver is hard to perform as soon as 
the sizes of source and recording device are not negligible compared to the wavelength. In 
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this paper, we propose that CBE can be observed without a standard source-receiver con- 
figuration. The opportunity to retrieve the CBE passively will be provided by analyzing the 
time-correlation of diffuse fields generated by a set of distant deterministic sources (exp. 1) or 
diffuse field generated by thermal phonons (exp. 2). 

The possibility of retrieving the impulse response (the Green function) of a complex 
medium by correlating records between two passive sensors has been suggested several times in 
the literature [f3]. Recently, Weaver and Lobkis [f4] proposed to cross correlate multiply scat- 
tered waves, a technique that intensively exploits the mesoscopic nature of diffuse wavefields, 
and has transpired to become a useful and accurate route to passive imaging applications 
ranging from ultrasound and ocean acoustics to seismology [15-f7]. In the following we apply 
this passive technique to observe CBE within an array of sources, the wave field being sensed 
by a distant receiver. By reciprocity, this configuration is analogous to an array of receivers 
and a distant source. 



Fig. 1 - First experimental setup: Correlation of actively generated diffuse fields. A series of 60 laser 
shots Si is performed on the top of a reverberant aluminum body, the ultrasonic wave field is sensed 
by a distant receiver R. CBE is observed within the Si's. 

The configuration of the first experiment (correlation of an actively generated diffuse field) 
is depicted in Fig. [TJ A piezoelectric transducer R is oil-coupled to the top surface of an 
irregularly shaped object whith dimensions '^90xl25xf25 mm. A hole is drilled through the 
block to break much of its symmetry and enhance the ray chaos. A Q-switched 60 mJ Nd-Yag 
infrared laser with pulse duration ~8 ns and a beam diameter ~f .5 mm excites elastic waves 
by means of ablation mechanism at position Si, resulting waves then reverberate thousands 
of times in the aluminum body. The receiver is far (several wavelengths) from the array of 
sources, its precise position does not play any role in our experiment. The absorption time 
Ta varies with frequency and is the order of tens of milliseconds, so the records are well above 
the noise level over the whole T=2f0 ms acquisition duration. The signal cpi{t) for the source 
i is then low-pass filtered, digitized and passed to the computer for processing. The laser 
source triggers the acquisition. The laser is mounted on a scanning motor of 1.5 mm step, the 
pitch-catch sequence is repeated over i = 1..60 positions. When the last impulse response is 
obtained, the laser is moved back to the initial position. To check the change of the medium 
due to temperature, one additional response (p'i{t) is acquired and compared to the initial 
signal (pi(t). Within the few minutes of the experiments, the thermal variations [f8] were 
found to have no effect on the results presented in this article. The time-averaged correlation 
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is processed for each couple of sources ij for times r ranging from to 60 ms: 

Q,(t) = / •f,{t)^^{t + T)dt. (1) 

Jo 

The configuration and processing of this first experiment are analogous to a one-channel 
time reversal experiment [f9], where the field is emitted in Si, sensed and time-reversed in 
TR, back propagates in the medium and is finally received in Sj. The analogy between 
time-reversal and correlations has been developed by Derode et al. [15] and holds as long as 
reciprocity is valid. A schematic view of the analogy is presented in Fig. [21 The link between 
time- reversal and CBE was also developed by De Rosny et al. [20]. In the open scattering 
medium they used, they found an energy enhancement of the time-reversed field if the source 
S and the receiver TR (the time-reversal device) are close together. As we will see, an en- 
hancement is also visible in a finite body around S even if R and S are far apart. 
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Fig. 2 - Three possible experimental configurations to observe CBE. (a) Conventional setup: Si is 
a source, the wave fields sensed in Rj after propagating in the medium, (b) Passive CBE: Si and 
Sj are two sources, the wave fields are sensed in 7? before correlation, (c) Time-reversal experiment: 
the field is emitted in Si and sensed in TR, then time-reversed and send back to the medium to be 
probed in Rj. 



The correlations are filtered with a / f =50% bandwidth, the intensity is averaged 
over all the ij couples available for a given distance r — Si — Sj, and additionally over 
times 30 ms< r <60 ms greater than the Heisenberg time of the cavity, which at 87 kHz is 
Th ~ 16 ms. The spatial enhancement of CBE is defined as the following ratio: 

CBEspaUai{r) = (^c^,lf-^x))' 

If the integration time T and the correlation time t are well above the Heisenberg time, 
then a decomposition of the wave field in uncorrclated modes Un is relevant and the ensemble 
averaged correlations reads [21]: 

(^') = {Y.u^{Sifu^{S,fu^{RY'^ (3) 

We assume that the receiver R is far from the array of sources Si-j, and that the modes are 
random gaussian variables with a spatial correlation of: {un{Si)un{Sj)) — Joi'^^r) where 
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Jo is the Bessel function of zero order. Therefore, after an expansion of Eq. [31 all the terms 
{un{R)un{Siorj)) — vanish. Each mode having the same mean-intensity (u^) , Eq.[3]rewrites: 

{CD = {u^f {{u'f + 2 {ur^SMS,)?) (4) 

Therefore the normalized CBE at a given frequency / averaged over the A/ frequency 
band takes the theoretical form: 

1 r^i+24fM,V/- (5) 




Fig. 3 - Spatial extent of the coherent backscattering effect for time r greater than the Heisenberg 
time. Experimental data are processed following Eq. (2) after a /=87 kHz 50% bandwidth filter. 
Theory comes from Eq. (4). 



Theoretical and experimental CBE are displayed in fig. [3] around f—87 kHz. The relevant 
velocity here is the Rayleigh wave speed c — 2.9 mm//is. The theoretical curve satisfactorily 
fits the experimental data, including the secondary lobes of the Bessel function. The passive 
setup allowed an optimal spatial resolution, which was only limited by the size of the source, 
and could be well bellow 1 mm if we chose to use a focusing optical lens. In an active source- 
sensor experiment, CBE is observed under the condition that sources and receivers have the 
same symmetry. This condition is naturally ensured here since we evaluate the impulse re- 
sponses between reproducible sources S. 

Let's now focus on the dynamics of the enhancement factor. In an active experiment the 
enhancement factor was found to vary with time t in the coda [6, 21] as 3 — b{T/TH), where 
b is the Fourier transform of the Dyson's two-level cluster function, b decreases from one at 
short times to zero at late times. Therefore the enhancement factor should increase from two 
for r < to three for r > Tff- 

Because correlations (r) are not exactly the impulse response between i and j but also 
have residual fluctuations resulting from incomplete ensemble or time averages, the enhance- 
ment factor and its dynamics are here slightly different from that of an active CBE. The 
intensity N{T) of these fluctuations of the diffuse field correlations has been studied [22], but 
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this theory does not hold for Si — Sj; a. theory at backscatter is needed. Fluctuations of the 
correlations of thermal noise do not depend on the distance between the receivers but only 
on the record length T. Therefore, instead of studying the dynamics of CBE with an actively 
generated diffuse field, we chose to realize another experiment based on correlation of thermal 
noise [23]. 



Fig. 4 - Second experiment: correlation of thermal noise. Two acoustic emission transducers sense 
the ultrasonic thermal noise at Ri and Rj. Waveforms are acquired after a 100 dB amplification, 
digitized and stored on a computer. 

The second experimental setup is depicted in Fig. 3) Because thermal noise is much weaker 
than the diffuse field obtained in the first experiment, we need to employ more sensitive 
acoustic emission transducers, noted Ri and Rj, along with two 100 dB pre-amplifiers (from 
Panametrics) . The distance Ri — Rj is much greater than one wavelength. The aluminum 
block has a size comparable to that of the first experiment. Thermal noise is recorded during a 
period of time T=400 s at a sampling frequency of 2.5 MHz; 6 series of acquisitions at different 
position pairs i ^ j are performed. The total amount of waveforms collected represents 
24 Gbytes of data. The useful thermal noise from the aluminum block represents only 40% of 
the record, but the additional electronic noise is white and (5-correlated, and therefore average 
down in the correlation. Correlations are then processed, and the dynamics (time-dependence) 
of the enhancement factor is evaluated through the following ratio, which compensates for the 
transducers' sensitivities [24]: 



We implicitly assume a spatial averaging over the six available positions, and a time- 
averaging in a 2 ms sliding window centered around r. Results are presented in Fig. [5] for 
two different frequencies (100 and 250 kHz) and two different record lengths (T = 10 s and 
T = 400 s). The Heisenberg time of the block around 100 kHz is Th — 18 ms, and is 
Th — 53 ms around 250 kHz. The fluctuations N{T), which are mainly determined by the 
amplifiers electronic noise, show a constant intensity in time r. They also linearly increase 
with the record length T whereas the intensity of the Green function in the correlations 
increases quadratically. Because of the large surface of contact of the transducers and the oil 
couplant, the absorption time describing the intensity decay is much smaller than in the first 
experiment: = 4.5 ms around 100 kHz and = 5 ms around 250 kHz. Therefore the mean 
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Fig. 5 - Dynamics of CBE obtained by correlating thermal ultrasonic noise. Experimental data are 
processed following Eq. [S] Theory without loss (absorption) is 3 — &(r) [21], and theory with loss is 
from Eq.[7][25]. 



intensity of the correlation is well described by Eoe^'^^'^" +N{T). At early time t, fluctuations 
are negligible: Eo/N{T) — 400 at r = 1 ms and for T = 400 s. As r increases, correlation 
intensity decreases, and fluctuations are finally dominating the correlations for times greater 
than 50 ms. We propose the following as a model for the passive CBE dynamics: 

rnF ,_. E,{i-h{T/TH))e-^/^^+N{T) 

^DrjdynarmcKT) - EqC-^/'^'- + N{T) 

This theory takes into account loss in the elastic body Ta and fluctuations in the correlations 
N{T). It is displayed in solid line in Fig. [5] and is found to well fit the experimental data 
(circles): for early t, CBE is increasing from two to three until fluctuations dominate and the 
enhancement factor drops to one. 

In this paper, we conducted two experiments on ultrasound propagating in an elastic body. 
In the first one, we retrieved the spatial shape of the coherent backscattering enhancement 
by correlating diffuse fields. This passive technique allows to observe CBE within an array of 
sources using only one or two distant receivers, and the reciprocity principle. Additionally, 
CBE shows a spatial resolution much greater than previous active experiments [6,8, 11]: the 
resolution here is no longer limited by the source-sensor distance but only by the distance 
between two consecutive sources. We also recalled the link between time-reversal, correlations 
and CBE. Godin [26] recently demonstrated that in a non-reciprocal medium, correlations still 
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yield the Green function. Therefore, as observed in an active experiment [20] , we suspect that 
breaking the reciprocity will destroy our paNsivc CBE. In the second experiment presented 
here, we correlated ultrasonic noise made of thermal phonons. The dynamics of CBE was 
studied and the main features were recovered: the increase of the enhancement factor from 
two to three is driven by the Heisenberg time. The effect of loss in the aluminum body and 
fluctuations in the correlations were also evaluated and proved to be a possible limitation. 
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